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Charge asymmetry in processes e''"e“ —> pA^ and e''"e“ —^ 7r"'"7r“7 is measured using 232 fb“^ 
of data collected with the BABAR detector at e^e“ center-of-mass energies near 10.58 GeV. An 
observable is introduced and shown to be very robust against detector asymmetries while keeping a 
large sensitivity to the physical charge asymmetry that results from the interference between initial 
and final state radiation. The asymmetry is determined as a function of the invariant mass of the 
hnal-state tracks from production threshold to a few GeV/c^. It is compared to the expectation from 
QED for e^e“ —>■ p~'y, and from theoretical models for e^e“ —>■ Tr+Tr'y. A clear interference 
pattern is observed in e''"e“ —>■ 7r^7r“7, particularly in the vicinity of the /2(1270) resonance. The 
inferred rate of lowest order FSR production is consistent with the QED expectation for e''"e“ -A 
p'^p~"f, and is negligibly small for e''"e“ -A n^'K~'y. 

PACS numbers: 13.40Em, 13.60.Hb, 13.66.Be, 13.66.Jn 


I. INTRODUCTION 

The radiative processes 

e'''e“ —?► Xj (1) 

have been extensively studied by several e“'"e“ exper¬ 
iments and the cross sections for e“''e“ —?> X have 
been measured using the initial state radiation (ISR) 
method W- At BaBAR , the cross sections have thus 
been determined in large energy ranges below the to¬ 
tal e+e“ center-of-mass (c.m.) energy ~ 10.58 GeV 
available at the SLAG PEP-II collider. The state X 
can be either fully described by Quantum Electrody¬ 
namics (QED) such as p'^p~, or any hadronic state with 
JPC = 1—. 

In reaction Q at lowest order (LO) the photon can be 
emitted from either the incoming electron or positron, 
or from the final state (final state radiation, or FSR). At 
BaBAR, the kinematic conditions are such that the pro¬ 
cess is dominated by ISR photons, which justifies the 
ISR method. The LO FSR contribution to the hadronic 
radiative process is neglected, as its theoretical estimates 
are well below the systematic uncertainties of the cross 
section measurement. This is due to the fact that the 
available e+e“ c.m. energy is far beyond the domain of 
the hadronic resonances that dominate the cross-section, 
so that hadronic form factors considerably reduce the 
probability that the photon is emitted from the final 
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state. However, the theoretical estimations are model- 
dependent, and it is thus important to have a direct 
experimental proof of the smallness of the FSR contri¬ 
bution to the hadronic cross sections when high preci¬ 
sion is at stake, as for the determination of the hadronic 
contribution to the g — 2 value of the muon . Because 
of the point-like nature of the muon, the FSR reduction 
does not occur for the e+e” —i pPp~^ process. The 
LO FSR contribution to the cross section is expected to 
vanish at threshold and to increase with the invariant 
mass of the muon pair (m^^). Still, the FSR fraction 
remains small for low di-muon mass (less than 1% for 
< 1 GeV/c^). For the e+e“ —)■ p^p~"f cross section 
measurement, a correction is applied for the LO FSR 
contribution as a function of which is so far de¬ 

termined by turning off FSR in the Monte Garlo (MG) 
generation. 

While it is not possible to distinguish ISR from FSR 
photons on an event-by-event basis, as the correspond¬ 
ing amplitudes are both present and interfere, a mea¬ 
surement of the interference provides a sensitive and 
quantitative determination of their relative strength. 
Measurement of the forward-backward asymmetry of the 
pions was first proposed in Ref. [5], as a test of the 
underlying model for final state radiation. In this pa¬ 
per, the ISR-FSR interference for e+e“ —>■ pPp~j and 
e+e“ -A- 7r+7r“7 is studied through the charge asym¬ 
metry of the production of these events at various de¬ 
cay plane angles. The comparison between the QED 
prediction and the measurement is done for the charge 
asymmetry in e“'"e“ -A pPp~j. Various FSR models 
are discussed for e+e“ —7r+7r“7, and the most realistic 
quark-FSR model is compared to the measurement of 
the charge asymmetry in that channel. 

This paper reports the first measurement of charge 
asymmetry in the e“'"e“ -A p^p~j process. For e“'"e“ -A 
7 r+ 7 r“ 7 , a preliminary measurement [7] of the forward- 
backward asymmetry has been reported at low energies 
{yAs ~ IGeV). No previous result exists at high ener¬ 
gies. 
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II. ISR-FSR INTERFERENCE AND CHARGE 
ASYMMETRY 

A. Charge asymmetry 

The Feynman diagrams for the LO ISR and LO FSR 
emission in the process e+e“ —> x'^x~^ (where x = fi oi 
tt), are illustrated in Fig. The total LO amplitude A4 
is the sum of the corresponding amplitudes AIisr and 
A^fsr, and the cross section for e+e“ —>■ x~^x~^ is 

a oc \M\^ = IAIisrP + l-^FSRp + 27?.e(A^isRA^FSR)- 

( 2 ) 

If the photon is emitted from the initial (final) state, 
the x~^x~ pair is produced with charge parity C = 
— 1 (+1), which implies that the interference term 
changes sign if one interchanges x"*" and x~. While the 
contribution of the interference term to the total cross 
section vanishes when one integrates over the kinematic 
variables of the final state, that term induces a signif¬ 
icant observable charge asymmetry in the differential 
cross section. 

Charge asymmetry is defined as 
^ \M\^ 

^ 27^e(AIISRAIFSR) /on 

IAIisrP + IAIfsrP ’ 

where x'^ -O- x~ means that x'^ and x~ are interchanged. 

Although it is not possible to reconstruct AIisr or 
AIfsr from the charge asymmetry and the cross section, 
as the relative phase between them remains unknown, 
information on the ratio IA^fer/AIisrI can be derived 
within the framework of specific models. 


B. Choice of kinematic variables 

Aside from an overall azimuthal rotation about the 
beam axis, the kinematic topology of the x+x “7 final 
state (where x = /x or tt) is described by four variables, 
which are the muon-pair (pion-pair) invariant mass mxx 
(or equivalently if*, the energy of the radiated photon 
in the e^e“ c.m.) and three angular variables. At a 
given nixx mass, the distribution of the three angular 
variables contains all the available information on the 
ISR/FSR amplitudes. 

At variance with the definition of forward-backward 
asymmetry used in Ref. [8], which refers to the polar 
angle of x~ with respect to the incoming electron in the 
e“'"e“ c.m. system (c.m.s.), this analysis introduces the 


set of angular variables illustrated in Fig. These are 
found to be more sensitive observables to measure the 
ISR-FSR interference: 


• 0* — polar angle of the radiated photon in the 
e+e“ c.m.s. (with respect to the e+e“ axis), 

• 0* — polar angle of x~ with respect to the photon 
axis in the x“'"x“ c.m.s.. 


• (jf — azimuthal angle of x~ with respect to the 
7 e“'"e“ plane in the x“'"x“ c.m.s. (or the e“'"e“ 
c.m.s.) 

Since x"*" O x~ interchange means reversal of the x~ 
direction to its opposite in the x+x” c.m.s. system, the 
charge asymmetry, for fixed nixx and 0*, is equal to 


^ g(0*,(^*)-a(7r-0*,7r-h0*) 

’ a{9*,(j)*) + a{Tr — 9*,Tr + (j)*) 


( 4 ) 


For the e+e“ —>■ process, the charge asym¬ 

metry as a function of cos0 * and </)*, studied with the 
AfkQed generator (see Sec. IVB), is shown in Fig. 
The FSR amplitude is dominant at |cos0*| ~ I, when 
one of the charged-particle tracks is very close to the ra¬ 
diated photon. However, Fig. shows that (/>* is a more 
sensitive variable to measure the ISR/FSR content over 
the full phase space, with sign reversal of the charge 
asymmetry. After integration over cos 9* and integra¬ 
tion over symmetrical cos 9* intervals, the distribution 
of the integrated charge asymmetry A(cos())*) suggests 
a simple linear dependence 


A(cos (/*) = Aq cos (j)*. 


( 5 ) 


From the expressions of the differential cross section de¬ 
tailed in the next section, it results that the slope Aq 
is an estimator of the ISR-FSR interference, sensitive to 
the ratio |A1fsr/AIisr| in each interval. Moreover, 
it will be shown in Sec. |V| that the measurement of Aq 
is barely affected by detector charge asymmetries. 


III. THEORETICAL PREDICTIONS FOR THE 
CHARGE ASYMMETRY 


A. QED prediction for the e“*“e —>■ 7 

process 


In the massless limit [^, the differential cross sec¬ 
tion of the QED e+e“ —)■ /x+/x “7 process, written as 
a function of the four kinematic variables defined above 


(Sec. II B I, implies that the differential charge asymme¬ 
try is proportional to cos cj)*: 
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FIG. 1: Feynman diagrams for e^e —^ x'^x -y {x = fj,, tt), where the photon is from lowest order initial state radiation (LO 
ISR, left) or lowest order final state radiation (LO FSR, right). 




FIG. 2: Definition of the angular variables describing the kinematic topology of the final states of the process e^e —>■ x^x 7 
{x = HyTr) at a given x^x~ invariant mass: (left) in the e"'"e“ c.m.s., (right) in the x'^x~ c.m.s. 



FIG. 3: Charge asymmetry at generator level in e"'"e“ —>■ simulation, as a function of cos 6* and <f)* for the same 

interval (6.5 < < 7.0GeV/c^) and various cos 6 '* ranges: (left) —1 < cos6* < —0.6, (middle) —0.6 < cos0* < —0.4, (right) 

—0.4 < cos 9* < 0. 








































































































































































































( 6 ) 


A.f._ 






1 = 


When the masses are taken into account, the effect 
from the electron/positron mass is found to be negligi¬ 
ble for radiated photons away from the beams. The 
effect from the muon mass is sizeable, especially at 
large when the radiated photon is close to one 
of the muons. Predictions for the charge asymmetry 
in the massive case are obtained by numerical integra¬ 
tion of several variants of the QED differential cross 
section HE]. The phase space considered in those 
calculations is limited to the experimental acceptance 
20° < 9* < 160°, and the results are shown in Fig. 
as a function of Predictions differ at the phys¬ 
ical threshold = 2m where only the charge 

asymmetry based on Ref. HH extrapolates to zero as 
expected, suggesting that the validity of formulae in 
Ref. H Uni does not extend to small At large 

mass > 3 GeV/c^), the prediction from Ref. [9] dif¬ 
fers from the others by up to a few percent. The formula 
of the differential LO cross section implemented in the 
AfkQed generator, which is used in this analysis for sim¬ 
ulation (see Sec. IVB I, is the one by Arbuzov et al. [m, 
which has the most reliable behavior over the full 
range. 


B. FSR models for the e“*“e —>■ 7r“*“7r 7 process 

As in the e+e“ —>■ process, ISR and FSR con¬ 

tribute to e+e“ —?> 7 r+ 7 r “7 (Fig.[^. However, the charge 
asymmetry is expected to be much smaller in the lat¬ 
ter process because the FSR contribution is strongly re¬ 
duced by the pion form factor at large y^. In addition, 
its estimate is model-dependent. 


1. FSR from point-like pions (model 1) 

In the FSR model shown in Fig. [^b), the photon is 
emitted from one of the final state pions, where the pion 


2y/s m^^ sin 0* sin 9* cos </>* 
s sin^ 9* + ml,, sin^ 9% 

MM 1 


is treated as a point-like particle. In this hypothesis, 
the FSR amplitude AIfsr is proportional to the pion 
form factor at the collision energy squared s, namely 
FV(s)- The ISR amplitude AIisr shown in Fig. [^a) is 
proportional to the pion form factor ^-t a reduced 

energy squared s' = s(l — 2E*/^/s). According to this 
FSR model, the charge asymmetry to be measured at 
BaBAR reflects the relative magnitude of the pion form 
factor at = 10.58 GeV and at low energy. It is con¬ 
sequently negligibly small, since Mr (s'), dominated by 
the p resonance in the s' = domain accessible to 
the experiment, is three orders of magnitude larger than 
|M.(10.58^ GeV^)| 0.01, as estimated from an extrap¬ 
olation of existing data mm using a 1 /s dependence. 
This model is studied with the PHOKHARA 4.0 [13] gen¬ 
erator, in which the FSR current has a point-like Lorentz 
structure, including a contact term, globally multiplied 
by the pion form factor. In this model, the Ao(m 7 r 7 r) 
distribution is expected to increase quadratically with 
mass on the p resonance, with a change of sign at the p 
mass 


Ao~2xl0-3(mL-m2), (7) 

with values well below the sensitivity of this analysis 
because of the large pion form factor suppression at 
10.58 GeV. 


2. FSR from quarks (model 2) 


In the a priori more realistic FSR model for e^e“ —>■ 
7 r+ 7 r “7 depicted in Fig.j^c), the FSR photon is emitted 
from the quarks, which subsequently hadronize into a 
pion pair El- The dominant ISR and FSR contribu¬ 
tions, and their interference, are written in terms of the 
variables defined in Sec. imi 


j ISR 

e+e~—>-7r+7r~7 




dml^ dcos9* dcos9* dcj)* 16Trs‘^ml^{s — m^^) 




.9 / 1 + cos^ dl „ ^ ^ 

(s^ -I- m^^) - . , „ ^ sin^ 9 4sm^„. cos^ 9 


sin^ 9* 


-\-2y/sm.,rr{s + m^„.)(tan0/) ^ sin20* cost/* — 2sm1^ sin^ 9* cos2(/*|. 


( 8 ) 


X 
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FIG. 4: (left) Charge asymmetry at (f)* = 0, 2 lo, as a function of obtained by numerical integration according to three 
different theoretical predictions (see text), with the condition 20° < 9* < 160° applied, (right) The difference between the 
prediction and the AfkQed LO value. Results labeled GW, BK, AF are obtained from references [9],[iQ],[n], respectively. 




(b)Final state radiation with pions 
treated as point-like particles (FSR 
model 1) 


7 



(c)Final state radiation at quark level 
(FSR model 2) 


FIG. 5: Feynman diagrams for e^e —>■ tt+tt 7 . 


where a and P are the QED fine structure constant and the pion velocity /? = a^I — respectively. The 

FSR contribution is 


e+e —>- 7 r^ 7 r 7 

dcos 9* dcos 9* d(p* 

and the interference term 


a^l3{s - ml^) 
647rs^ 


{l + cos^9;)\V{mi^,9*)\\ 


+ e “ —f TT + TT “ 7 


dm^^ dcos 9* dcos 9* dcj)* 


-Re{F:{ml^)V{ml^,9*)} 


I r- r/i 2 2 . 2 /,*iSin0* cosd)* 

X cos 9^ cos 9 -F [(I-F cos^ 0.^)5-F sin^ — . - — 

z sm 


where 


r= “S' 


dz 


2z- I 

^(l - z) 


®|( 2 ,wL>cosr) {q = u,d), 


(9) 


(lO) 


(ll) 


and $+(z,TOj^,cosd*) is the C-even part of the 2-pion time-like form factor is taken from a fit to 

generalized distribution amplitudes (GDA). The pion 
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BaBAR data [S] with a vector dominance model. 

So far, there is no implementation of this model in 
an MC generator to describe the ISR-FSR interference 
in the e“'"e“ -A- 7r+7r“7 process. In order to predict 


the charge asymmetry numerically, we take the following 
GDA model, which is a modified version of the model 
found in Ref. [15]: 


^Ziz,m^^,cose*) 


10 z(l - z)( 2 z-l) 



where cq and C 2 are the magnitudes of the S-wave and 
D-wave contributions, respectively. As the scalar sector 
is known to involve wide resonances, the S-wave contri¬ 
bution is approximated by a constant amplitude with 
a mass-dependent phase do(™ 7 r'n-) taken from pion-pion 
phase-shift analyses m in the region below 1.6 GeV/c^. 
This model incorporates the rapid phase variation across 
the /o(980) resonance. Using cg = —0.5 [TSI yields an 
Aq value of about —1% near the p resonance and nearly 
flat with mass. For the D-wave tensor contribution, we 
use a Breit-Wigner form (BW) for the /2(1270) reso¬ 
nance in order to take properly into account the mass 
dependence of the amplitude, the phase variation being 
given by the BW form in agreement with the measured 
52 {‘niT^Tr) values [l6|. The angular dependence in the 
TTTT center-of-mass is given by the Legendre polynomial 
P 2 {cos 6 *), which assumes the dominance of helicity 0 
for /2(1270) production. 


C. Other sources of charge asymmetry 

Next-to-leading order (NLO) corrections including ad¬ 
ditional photons (soft and hard) and loops are expected 
to affect the lowest-order (LO) predictions for the charge 
asymmetry. For the pT^ process these corrections 
have been computed recently m and implemented in 
the PHOKHARA 9.0 generator [TH|. As discussed in 
Sec. |VIG| the effects are found to be small, at the per¬ 
cent level for the experimental conditions of the present 
analysis, and to be well accounted for by the simpler 
structure function approach implemented in AfkQed. 
No exact NLO calculation is available for the 7r+7r“7 
process. In this case, since the LO charge asymme¬ 
try is expected to be small because the FSR ampli¬ 
tude is suppressed, NLO corrections could play a rela¬ 
tively more important role. The soft and virtual photon 
contributions to the Born process e^e“ —>■ tt+tt” are 
known mnoi to generate an asymmetry of the pion 
production, with asymmetry values at the percent level 
at a TTTT mass of IGeV/c^. However, it is unclear if the 
above result can be used in the conditions of the present 
process e“'"e“ —>■ 7r“''7r“7, where one of the incoming 
electrons is highly off-shell after emission of a hard ISR 




r 


+c2/3"BW(m^^)P2(cosr) 


( 12 ) 


photon. Furthermore, such an asymmetry would vanish 
because of the symmetrical integration in cosd*. NLO 
corrections as implemented in AfkQed have indeed no 
effect on the charge asymmetry. No correction on the 
measured charge asymmetry is therefore applied for the 
7 r+7r“7 process. 

Another potential source of charge asymmetry comes 
from Z exchange. This contribution is strongly sup¬ 
pressed by the Z propagator, especially for the ISR di¬ 
agrams where ~ 10“^. Therefore one expects 

this effect to be negligible for 7r+7r“7. The contribution 
is larger for the FSR diagrams for p'^p~"/ since here the 
relevant ratio is s/M"^ = 1.4%. The contribution of Z 
exchange is studied with the KKMC generator |3D|. As 
the effect is at the level of a few 


reported in Sec. |VI G 
per mille. 


IV. EXPERIMENTAL ANALYSIS 
A. The BABAR detector and data samples 

The analysis is based on 232 fb~^ of data |2T| collected 
with the BaBAR detector at the SLAG National Accelera¬ 
tor Laboratory at the PEP-II asymmetric-energy e“''e“ 
collider operated at the T(4S) resonance. About 10% 
of the data was collected 40MeV below the resonance. 
The BaBAR detector is described in detail elsewhere [52] . 
Gharged-particle tracks are measured with a five-layer 
double-sided silicon vertex tracker (SVT) together with 
a 40-layer drift chamber (DGH), both inside a 1.5 T 
superconducting solenoid. Photons are assumed to orig¬ 
inate from the primary vertex defined by the charged- 
particle tracks of the event, and their energy and posi¬ 
tion are measured in a CsI(Tl) electromagnetic calorime¬ 
ter (EMC). Charged-particle identification (PID) uses 
the ionization energy loss dE/dx in the SVT and DGH, 
the Cherenkov radiation detected in a ring-imaging de¬ 
vice (DIRC), the shower energy deposit (Ecai) in the 
EMC, and the shower shape in the instrumented flux 
return (IFR) of the magnet. The IFR system is made of 
modules of resistive plate chambers (RPC) interspaced 
with iron slabs, arranged in a layout with a barrel and 
two endcaps. Collision events are recorded and recon- 
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structed if they pass three levels of trigger (hardware, 
online software, and offline filter), each using comple¬ 
mentary information from the sub-detectors. 


B. Monte Carlo generators and simulation 

Signal and background processes e+e“ —>■ are 

simulated with the AfkQed event generator, which is 
based on QED for e+e“ —>■ and Ref. [53] for 

hadronic production. LO ISR and FSR emission is sim¬ 
ulated for e+e“ —)■ /i+/r“ 7 , while LO FSR is neglected 
for hadronic processes. The main photon (hereafter 
called ‘ISR’ photon) is emitted within the angular range 
20° < 9* < 160° in the e+e“ c.m. system, bracket¬ 
ing the photon detection range with a margin for res¬ 
olution . Additional ISR photons are generated with 
the structure function method [53], and additional FSR 
photons with the PHOTOS [25] program. Additional 
ISR photons are emitted along the e'^ or e~ beam par¬ 
ticle direction. A minimum mass tox 7 isr ^ 8 GeV/c^ 
is imposed at generation, which puts an upper bound 
on the additional ISR photon energy. Samples corre¬ 
sponding to 5 to 10 times the data are generated for the 
signal e^e“ —>■ and e^e“ —>■ 7 r''' 7 r “7 channels, 

as well as large samples of backgrounds from the other 
two-prong and multi-hadron ISR processes. Background 
processes t 99 (9 = u,d,s,c) are generated with 

the JETSET [55] generator, and e+e” —>■ r+r“ with the 
KORALB [27] program. The response of the BaBAR de¬ 
tector is simulated using the GEANT4 [55] package. 


C. Event selection 

Event selection follows the same procedure as the se¬ 
lection of two-charged particle ISR events used for cross 
section measurements [ 6 ]. It requires a photon with en¬ 
ergy E* > 3GeV in the e“'"e“ c.m. and laboratory polar 
angle with respect to the e~ beam in the range [0.35- 
2.4] rad, and exactly two tracks of opposite charge, each 
with momentum p > IGeV/c and within the angular 
range [0.40-2.45] rad. If more than one photon is de¬ 
tected, the candidate with the highest E* is taken to be 
the ‘ISR’ photon. To ensure a rough momentum balance 
at an early stage of the selection, the ‘ISR’ photon is re¬ 
quired to lie within 0.3 rad of the missing momentum 
of the charged particles (or of the tracks plus the other 
photons). The tracks are required to have at least 15 


hits in the DCH, to originate within 5 mm of the colli¬ 
sion axis and within 6 cm from the beam spot along the 
beam direction, and to extrapolate to the DIRC and 
IFR active areas in order to exclude low-efficiency re¬ 
gions. Both tracks are required to be identified either 
as muons or as pions. To suppress the background to 
7 r+ 7 r “7 at threshold due to the e“''e“ —>■ 77 process fol¬ 
lowed by a photon conversion and misidentification of 
both electrons as pions, it is further required that the 
distance in transverse plane Vxy between the vertex of 
the two tracks and the beam collision point be less than 
0.5 cm for m,r 7 i- < 0.5GeV/c^. Electron background to 
fi~"f is negligible over the full mass range. 

In order to suppress multi-hadron ISR events and re¬ 
duce higher order radiative processes, the selected two- 
prong candidates are subjected to a one-constraint kine¬ 
matic fit to the e+e“ —>■ hypothesis {x = /i,7r), 

in which only the two good charged-particle tracks are 
taken as input and the corresponding missing mass is 
constrained to the null photon mass. The value of 
the kinematic fit is required to be less than 15. 


D. Charge asymmetry calculation 


For a complete topology of the final states, the az¬ 
imuth cj)* defined in Sec. |II B| should cover the 27r range. 
However, the event sample with x~ azimuth G [0, tt] 
is complementary to the sample with x'^ azimuth S 
[ 0 , 7 r], since (/>!]_= tt -1- </)! (mod27r) in every event. This 
allows to restrict to the range [ 0 , tt] with no loss of 
phase space. After integrating over 9* and 0*, the total 
event sample in a fixed vixx interval subdivides into two 
subsamples: one with (jf_ G [ 0 , 7 r] {N-), the other with 
e [0,^] (Ar+). 

We obtain separately the distributions in coscj)* of 
the 77^*°® samples in data, namely Ai°’°®(cos 0 *) with 
(f)- = (!>* and N^^{cos(j)*) with = (j)* . Distribu¬ 
tions of background events N^*^{cos </>*) are determined 
separately for each subsample, as described below. Like¬ 
wise, efficiencies are split into e±{cos cj)*) and computed 
using the full simulation of e+e” —)■ x'^x~x (x = /i,7r) 
events, with corrections for the differences between data 


and simulation (see Sec. VIG and VIIC). In a given 


W/i/i (jnTTTr) interval, the asymmetry at a given cos(p* 
is derived from the difference between the N-(coscj)*) 
and N+{cos(j)*) yields, corrected for efficiency, and is 
obtained from the following expression 


^ N_(coscl)*)/e-{cos(j)*) - N+{cos (j)*)/€+(coscj)*) 

N-(cos cj)*)/e-(coscj)*) -|- N+(coscj)*)/e+(coscj)*) ’ 


(13) 


where N±(coscj)*) = (coscj)*) — (coscj)*). 


Note that distributions of N±(cos cj)*) can be obtained 
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in each {m^x, 9*, cos 9*) cell of the phase space, and the 
asymmetry defined by Eq. (131 can be calculated. How¬ 
ever, this one-dimensional quantity H(cos^*) is a valid 
definition of charge asymmetry only when the variable 
cos 9* is integrated within a symmetric range. This is a 
consequence of the fact that the x'^ O x~ interchange 
means both (jf —)■ tt -|- 0* and cosd* —>■ — cos0*, and 
therefore a non-null value of H(cos^*) in an arbitrary 
cos 9* interval is not an intrinsic signature of ISR-FSR 
interference. 


It is worthwhile to note that the kinematic selection 
in itself is charge-symmetric. Hence the observed bias 
on the measured raw charge asymmetry is a cross effect 
of physical charge-asymmetric kinematics and charge- 
symmetric detector acceptance. It does vanish for a null 
physical charge asymmetry. As checked with a 
simulated sample produced by ISR only, no fake charge 
asymmetry emerges from the kinematic selection. 


2. Software trigger and tracking 


E. Event samples and backgrounds after selection 

The cos(j>* distributions AI2'^®(cos </)*) for e“'"e“ —>■ 
obtained in data after the overall event selec¬ 
tion are shown in Fig. for /x/r mass intervals rang¬ 
ing from threshold to 7 GeV/c^. The event distributions 
N^‘^{cos(j)*) obtained for e+e“ —>■ 7r+7r“7 in data are 
shown in Fig. Ih in O.lGeV/c^ mass intervals ranging 
from 0.3GeV/c^to 1.8GeV/c^. 

The backgrounds remaining after selection are es¬ 
timated using the full simulation, normalized to the 
data luminosity, of the non-signal two-prong ISR events, 
multi-hadron events produced through ISR, e+e” —)■ qq 
events, and events. The expected contamination 

for e+e” —)■ a function of cos cj)* in typical 

rrififj, intervals is shown in Fig. where the total error is 
the quadratic sum of the statistical error and 10% sys¬ 
tematic uncertainty on normalization [B]. Likewise, the 
estimated backgrounds for e“''e“ —>■ 7r+7r“7 in typical 
intervals are shown in Fig. 

V. ACCEPTANCE AND DETECTOR 
EFFICIENCY EFFECTS ON THE CHARGE 
ASYMMETRY 


Biases on the raw charge asymmetry measurement 
originate from the software trigger and the track re¬ 
construction. They are observed in the low mass re¬ 
gion, as illustrated in Fig. and vanish at high mass 
> l.SGeV/c^). 

The common origin of the mass-dependent trigger 
and tracking inefficiencies is geometrical and has been 
thoroughly studied for the tt+tt” cross section measure¬ 
ment [6]. Converging trajectories in the DCH of op¬ 
positely deflected tracks emitted in close-by directions 
confuses the track reconstruction and causes both the 
software trigger and the final tracking inefficiencies. In 
the charge-conjugate configuration, in which the positive 
and negative tracks are interchanged, tracks diverge in 
the magnetic field and are well separated in the trans¬ 
verse plane, although with the same absolute azimuthal 
opening angle. The efficiencies are consequently charge- 
asymmetric, sharply reduced for overlapping tracks in 
low regions, at Aej) values close to zero but always 
positive, where A(j) is the signed angular difference be¬ 
tween the azimuths of the positive and negative tracks 

A0 =((/)+-(^_) e [-7r,7r]. (14) 


The charge asymmetry measurement is affected by the 
event reconstruction and selection. These experimental 
effects are investigated using the full simulation of signal 
events through changes of the raw charge asymmetry, 
defined as = {N_ — N^)/{N_ + Ab|_), which are 
observed after each selection step as a function of cos (f* ■ 

A. Study of the effects with the muon simulation 

1. Kinematic acceptance 

The kinematic acceptance includes the angular ac¬ 
ceptance for the primary photon and the two charged- 
particle tracks, and the momentum restriction (p > 
1 GeV/c) applied to charged-particle tracks. Each kine¬ 
matic selection is found to modify the slope of the raw 
charge asymmetry significantly, though the total effect 
on the slope from the kinematic requirements altogether 
turns out to be small due to accidental cancelations. 


3. ‘ISR ’ photon reconstruction 

The event selection requires that an TSR’ photon with 
F* > 3 GeV be measured in the EMG. The raw charge 
asymmetries for the fully simulated events with and 
without the requirement of the ‘ISR’ photon being re¬ 
constructed are shown in Fig. El where the events are 
already required to be within the kinematic acceptance. 
Effects are observed in high regions. 

The origin of a charge-asymmetric photon reconstruc¬ 
tion inefficiency is again geometrical. In case one of 
the charged-particle tracks and the ‘ISR’ photon overlap 
in the EMC, the shower produced by the ‘ISR’ photon 
is mistakenly associated to the charged-particle track, 
and the ‘ISR’ photon is lost. In the charge-conjugate 
configuration, no overlap occurs because of the oppo¬ 
site deflection of the charged-particle track in the mag¬ 
netic field. The overlap happens at Acj) = —tt -|- e 
where e is a small positive quantity, and as a con¬ 
sequence, the ‘ISR’ photon reconstruction efficiency is 
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charge-asymmetric, strongly reduced around Ac/) ~ —tt. 
As the overlap of one charged-particle track and the 
‘ISR’ photon occurs preferentially at high mass, due to 
phase space, the corresponding effects are only observed 
at m^fj, > 3.5GeV/c^. 


4- Muon identification 

The charge asymmetry measurement is also affected 
by muon identification, in the low and high re¬ 

gions, as shown in Fig. |12[ Charge-asymmetric ineffi¬ 
ciency of muon identification results again from event 
topologies. The first cause, which affects low re¬ 
gions, is the overlap of the two charged-particle tracks 
at the IFR, which confuses the muon identification al¬ 
gorithm. The second cause is the partial overlap of one 
charged-particle track and the ‘ISR’ photon at the EMC, 
which makes the track look unlike a muon. The latter 
effect is more pronounced at high mass, when a muon 
and the ‘ISR’ photon are emitted in close-by directions. 
The efficiency of muon identihcation as a function of At/) 
exhibits a sharp dip at positive At/) at low mass, and at 
At/) ~ —TT at high mass. 


5. Summary of the acceptance and detector efficiency 
effects in the process 


ducing an asymmetry are nearly symmetric in cos (jf . 
The Ao observable introduced in Eq. ([^ is thus ex¬ 
pected to be robust against such effects. Overall, the 
slope of the asymmetry is barely affected by detector 
inefficiencies and event selection. In the yC 1 Pro¬ 
cess, the maximum effects, of a few 10“^, take place 
around ^ 2 — 4CeV/c^. In the low mass region 
[0.5-1.0] CeV/c^, the effect from the overall selection is 
at the level of a few 10“^. 


B. Study of the effects with the pion simulation 


The acceptance and detector effects are also studied 
with the simulated e^e“ —>■ 7r"''7r“7 events. The overall 
effect around the p resonance G [0.4, 1.2] CeV/c^) 
is (0.30 ± 0.07) X 10“^ in average. 

7r+7r“7 


As the charge asymmetry is null for e“'"e“ 

MC events, generated with no LO FSR, the accep¬ 
tance effects on the slope of the charge asymmetry are 
quite small for any selection requirement, including the 
kinematic ones. This is in contrast with the pApT^ 
case, where the individual kinematic requirements in¬ 
duce large effects. However this conclusion holds only if 
the charge asymmetry in the data is actually null. If a 
sizeable asymmetry is measured, the bias introduced by 
the cross effect of acceptance and asymmetry has to be 
evaluated and corrected (Sec. VIIB). 


The overall efficiencies e± needed to correct the N± 
event yields entering the charge asymmetry measure¬ 
ment (Eq. @) are the overall result of the acceptance- 
induced and detector asymmetries discussed above. 
They are determined using the full simulation, sepa¬ 
rately for the N± samples. 

As previously discussed, the detector inefficiencies are 
mostly caused by the spatial overlap of trajectories oc- 
curing in the detector: 2-track overlaps in the DCH and 
the IFR, respectively for Ai/i = O.I ± O.I and 0.5 ± 0.2 
and affecting masses below 2CeV/c^, and the photon- 
muon overlap in the EMC, for A^ > —tt and affecting 
masses above 4CeV/c^. These various overlap effects 
contribute very asymmetrically to the two N± samples, 
due to a complete correlation between the cost/)^- and 
A(/) variables. This is demonstrated in Fig. which 
shows that the N+ (V_) sample corresponds to Ai/) > 0 
(A(/) < 0). Since the 2-track overlaps occur for At/) > 0, 
one expects e_ > e+. For the photon-muon overlap with 
nearly opposite tracks the situation is reversed. 

As a consequence, as summarized in Fig. [Td] (Sec. [VI]), 
the acceptance and detector inefficiencies induce a 
change AA in the observed charge asymmetry magni¬ 
tude and also distort the linear dependence on cos/fi*. 
The dominant effects are from geometric acceptance, 
‘ISR’ photon reconstruction and the track momentum 
requirement p > I CeV/c. 

However, although the detector is not completely 
charge symmetric, Fig. shows that the effects pro- 


VI. RESULTS ON THE CHARGE 
ASYMMETRY IN THE e+g- ^ p+p~~f 
PROCESS 

The measured raw charge asymmetry after the com¬ 
plete event selection for the data is obtained as a func¬ 
tion of cos (/)* in various intervals, and shown in 
Fig. It is consistent to the first order with the 
full simulation of e+e“ —>■ p'^p~j events, except in 
the mass interval populated by the J/tf resonance (3.0- 
3.5CeV/c^). This is expected as J/if production is not 
considered in AfkQed. 

The physical charge asymmetry for e+e“ —>■ p^p~^ 
is obtained from the measured cos (/)* distributions after 
background subtraction and efficiency correction. The 
background dependence on cos (/)* is estimated with the 
simulation for each of the samples (j)’^ = (j)* G [0,7r], 
as explained above. Similarly, the overall efficiency 
e± is obtained with fully simulated e+e” —)■ p'^p~^ 
events, and corrected for data/MC differences in detec¬ 
tor response. The efficiency differences between data 
and simulation have been studied extensively for the 
cross section measurements of e+e” — p'^p~'y and 
e+e“ —>■ 7r+7r“7 [6]. They are parameterized as a func¬ 
tion of the azimuthal opening angle A(j) between the two 
muons, and projected onto the coscf* variable by sam¬ 
pling with MC. 

The charge asymmetry distributions for e^e“ —>■ 
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FIG. 11: Raw charge asymmetry as a function of cos0* in selected intervals, for e"''e“ ^ MC events with (o) and 

without (•) the requirement of the reconstruction of the ‘ISR’ photon, where the events are already required to be within the 
kinematic acceptance. 
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FIG. 12: Raw charge asymmetry as a function of coscfi* in selected intervals, for e''"e“ —> MC events with (o) and 

without (•) the two-muon identification, where the events are already required to be within the kinematic acceptance and have 
two good charged-particle tracks reconstructed. 


7 data after background subtraction and efficiency 
correction, as well as the charge asymmetry for e“'"e“ ^ 


/i^/r 7 MC at generation level, are shown in Fig. 15 


While the cos (/>* dependence of the measured raw charge 
asymmetry is not linear, the corrected data distributions 
are quite consistent with the MC distributions at gener¬ 
ation level. The slopes of charge asymmetry in various 
mass intervals are obtained by fitting the background- 
subtracted efficiency-corrected charge asymmetry distri¬ 
butions to ^0 cos (p*. 


A. Test of the charge asymmetry with 
J/tp —t events in data 

Since e+e“ —)■ "fJ/ip with on-shell J/ip —)■ is a 

pure ISR process at T(4S') energies, the J/ip —>■ 
sample in the data provides a test of fake asymmetries 
that could arise in the analysis. 


To overcome the limited statistics of the J/ip ^ 
sample, a loosened event selection is applied, with the 
muon identification requirement removed, which pro¬ 
vides a gain in statistics by a factor of about 4, with 
no significant increase of the hadronic background. The 
spectrum shows a clear J/ip peak, over a linear 
QED background. Defining and the re¬ 

spective slopes of the charge asymmetry for J/pj and 
underlying QED events, the slope Aq measured in the 
vicinity of the J/ip resonance is the average 


^j/tp + Nqed 


(15) 


where and Vqed are the yields from J/ip and 

QED, respectively. The quantities IVqed and Nj/^ are 
obtained by fitting the mass spectrum with a sum of a 
linear QED component and a Gaussian J/ip signal, with 
fixed width equal to the mass resolution at the J/ip and 
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0.5-lGeV 





FIG. 13: (left) Distribution of cos 01 vs. A0 for the N- sample (01 € [0,7r]) in the (0.5-1.0) GeV/c^ mass interval for the 
MG; (right) the same for the N+ sample (0+ € [0,7r]) 


centered at the nominal J/i!) mass. The slope 
is obtained by htting the charge asymmetry in the J/'0 
sidebands. The measured slope of the charge asymmetry 
as a function of is shown in Fig. 16 the expected 
behaviour is clearly observed, with a smooth variation 
from the QED continuum with a large negative value 
and a sharp peak approaching a null slope on the J/ip 
resonance. The specific slope is obtained as a func¬ 


tion of in Fig. [1^ according to Eq. (15): its value 


is stable across the j/'0 peak and a fit to a constant 
between 3.07 and 3.12GeV/c^ yields 


= (0.3 ± 1.6) X 10“ 


(16) 


which is consistent with zero, as expected from the ISR- 
only J/ip production. 


B. Comparison to QED 


10“^ at mass less than IGeV/c^ and above 5GeV/c^, a 
small but significant discrepancy shows up for intermedi¬ 
ate mass, reaching ~ 3 x 10“^ between 1.5 and 4 GeV/c^. 
Investigations of systematic uncertainties, both at the 
experimental and theoretical levels, are reported in the 
next section. 

Due to the asymmetry of the beam energies at PEP- 
II, independent charge asymmetry measurements in two 
different kinematic regimes are provided by splitting the 
data into a forward (cos 0* > 0) sample and a backward 
(cos0* < 0) sample. The full analysis, including back¬ 
ground subtraction and efficiency correction, is redone 
on each sample separately. The results are shown in 
Fig. A significant discrepancy between data and 
AfkQed is observed in the forward region, in the 1.5- 
4GeV/c^ mass region, while in the backward hemisphere 
data and AfkQed are consistent. The differences are 
quantihed in Table [ij No significant forward-backward 
difference is expected from the generator. 


The final slope Aq as a function of measured on 
the data is shown in Fig. |17[ together with the asymme¬ 
try at the MG generation level, and the difference be¬ 
tween them. The mass interval containing the pure-ISR 
contribution from the J/tp (e+e“ — » 7 isr J/V’i J/ip 




discussed in detail in Sec. VIA 


is excluded. 


The absolute difference between data and MG AAq = 
^data _ ^MC jg g percent level (0-3%). 

The measured slope of charge asymmetry is negative 
throughout the mass range under study, and its mag¬ 
nitude increases with mass, reaching values as large as 
—0.7 at 5 GeV/c^, in agreement with the trend predicted 
by QED. However, while data and LO QED agree within 


C. Systematic uncertainties 

1. Experimental systematic effeets 


The primary sources of systematic uncertainty are the 
background estimation, data/MG differences in detector 
response, and differences between the physical charge 
asymmetry in the data and in the generated MG events. 

The difference AAq between the results with and 
without background subtraction is found to be well be¬ 
low 10“^ except in the 0.5-1.0 GeV/c^ range, where it 
reaches 2x 10“^ because of the larger p background with 
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FIG. 14: Raw charge asymmetry as a function of cosi^* for fj,^fi 7 events in data (•) and MC (o), in various 
after the complete event selection. 
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FIG. 15: Charge asymmetry in e''"e —>■ 7 data before (•) and after (o) background subtraction and efficiency corrections, 

and for MC (A) at generation level. 
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FIG. 16: (left) The measured charge asymmetry as a function of mass near the J/iji resonance, (right) The derived charge 
asymmetry for from e''"e“ —>■ 'yJ/ip 




FIG. 17: The slope Aq of charge asymmetry as a function of in e^e“ —>■ data and in MG at generation level 

(left), and the absolute difference between them (right). The J/tp mass region (3.0-3.2GeV/c^) is excluded. Only statistical 
uncertainties are shown. 


TABLE I: The difference (in 10~^ units) between the measured Ao and the AfkQed prediction for the two mass intervals 1.5-4 
and 4-7 GeV/c^ in different cos 9* regions. The last line gives the difference between the two regions cos 9^ > 0 and cos0* < 0. 
Statistical uncertainties only. 


(10-2) 

1.5 < < 4GeV/c2 

4 < < 7 GeV/c^ 

all cos 9* 

cos 9* > 0 
cos 9* < 0 

difference 

2.65 ±0.38 

3.61 ±0.50 

1.08 ±0.60 

2.50 ±0.78 

0.86 ±0.22 

0.76 ±0.30 

0.82 ±0.31 

-0.05 ±0.42 
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Ag(data)-Ag(MCTmth) 



FIG. 18: The difference between the measured asymmetry 
slope in the process and the AfkQed prediction, as 

a function of excluding the J/'i/) 3.0-3.2 GeV/c^ region. 
Forward {cos 9* > 0) and backward (cos0^ < 0) hemispheres 
are analysed separately. Statistical uncertainties only. 


two pions misidentified as muons. Since the background 
level is known with better than 10% accuracy [H], the 
corresponding systematic uncertainty on the asymme¬ 
try slope is at most 2x 10““^ throughout the studied mass 
range. 

The trigger, tracking and ^-ID induce charge- 
asymmetric data/MC corrections, as overlap effects are 
not perfectly reproduced by simulation. However, the 
data/MC corrections have small effects on the charge 
asymmetry slope, at most 3x10“^. Since the correc¬ 
tions have been measured with a precision of 10% or 
better [S], the corresponding systematic uncertainty on 
the asymmetry slope is less than 3x10“^. 

As explained in Sec. |V A 1[ the effects from the kine¬ 
matic acceptance on the measured slope of the charge 
asymmetry depend on the physical charge asymmetry 
itself. The possible bias on the acceptance correction, 
induced by the physical charge asymmetry in the gen¬ 
erator inaccurately reproducing the data, is studied us¬ 
ing a sample of re-weighted e“''e“ —?► MC events 

where weights are adjusted to yield the same asymme¬ 
try as measured in the data in each {m^^, 0/, 0*, cos(/*) 
phase-space cell. The expected bias in the measurement 
from a difference of charge asymmetry between data and 
MC is found to be less than 5x10“^, which is taken as 
a systematic uncertainty on the Aq measurement. 


2. Effects from imperfect simulation 


Since a simple linear fit A(cos(/)*) = Aocoscj)* might 
be questionable, we perform an alternate two-parameter 
fit on the charge asymmetry after efficiency corrections 
A(coS(/*) = AqCOSi/* -I- Bq. The Bq values obtained 
in data are a few 10“^ at most, while the asymmetry 
slopes Aq deviate from the final values, which use the 
one-parameter fit, by less than 10““^. 

To investigate whether the observed discrepancy re¬ 
sults from the efficiency corrections, we study the differ¬ 
ence AA"'®''"(cos (/*) between the raw asymmetries ob¬ 
served in data and MC after full event selection. Al¬ 
though the raw asymmetry itself is not linear, especially 
at low mass (Fig. [I4|), the difference AA’'™(cos0*) in 
each mass interval is observed to be linear with coscf*. 
In particular, there are no edge effects in the vicinity of 
I cost/*I ~ 1, which could have resulted from different 
resolutions in data and MC. The results of the fits are 


shown on Fig. 19 (black points). The values of the slope 
of AA''™ are insensitive whether the linear fit is a one- 
parameter or a two-parameter fit. The data-MC discrep¬ 
ancy in the 1.5-4 GeV/c^ mass region is already observed 
at the raw level, which excludes efficiency or resolution 
bias. The Bq values returned by the two-parameter fit 
over IS.A'^™ depart from zero by up to 2.5x10“^ at low 
mass, as expected from imperfect detector simulation. 
When data/MC corrections of detector efficiencies are 
applied to the simulated raw data, the results are shown 
on Fig. 19 (blue triangles). The Bq values are reduced 
to a few 10“^, while the slope of A A''®’'" is not changed 
by more than 3x10“^. 

When performed on the forward and backward sam¬ 
ples independently, the study at the raw data level con¬ 
firms that the data-MC discrepancy in the 1.5-4 GeV/c^ 
mass interval is confined to the forward region, where 
the slope of AA''™ is significantly non-null by ~ Sdstat 
while it is consistent with zero within 2crstat over the full 
mass range for the cos0/ < 0 sample. In contrast, the 
fitted Bq values are consistent with each other in the two 
samples, except at very low mass (m^^ < IGeV/c^). 

Comparisons of the data and MC distributions of 
event variables entering the asymmetry analysis are per¬ 
formed, in particular near the acceptance boundaries. 
A sizeable departure is observed in the high E^, low 
9^ region, in the forward cos 9* > 0 hemisphere. How¬ 
ever, the asymmetry measurement is found to be insen¬ 
sitive to this discrepancy at very forward photon angles. 
To investigate whether different resolutions in data and 
MC might bias the efficiency corrections and the event 
assignment to the N+ or A_ samples, the analysis is 
fully redone with tighter acceptance requirements. The 
change of asymmetry slope is small, (4 ± 3) x 10“^. 
Conservatively, a systematic uncertainty of 7x10“^ is 
assigned to account for imperfect simulation near the 
edges of the selected phase space. 

The studies above first show that, although accep¬ 
tance and detector inefficiency effects are important. 
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FIG. 19: (left) Slope of the difference {cos (j>*) for between the raw asymmetries observed in data and in MC 

before (black points) and after (blue triangles) the data/MC corrections, as a function of (right) Constant term of the 

two-parameter linear fits of (cos <()*). Statistical uncertainties only. 


they are well accounted for in the simulation. Data/MC 
corrections are found to significantly reduce the symmet¬ 
ric component of the asymmetry {Bq), but most impor¬ 
tantly, the studies demonstrate that the measurement of 
the asymmetry slope is robust against uncertainties in 
the efficiency corrections. 

As a global test to differentiate between an uncor¬ 
rected experimental bias and a true deviation from the 
QED prediction, the difference between the measured 
asymmetry and the theoretical one, as implemented in 
AfkQed, is studied as a function of cos (j)* for the events 
in the mass interval 1.5-4 GeV/c^ where the deviation is 
the largest (excluding the J/V’ 3.0-3.2GeV/c^ region). 
As shown in Fig. a linear dependence is indeed ob¬ 
served. This supports the assertion that the deviation 
we observe does not originate from a detector effect un¬ 
accounted for in the simulation. 


3. Theoretical systematic effects 

The AfkQed event generator only includes the LO 
QED interference between ISR and FSR amplitudes. 
Additional photons generated independently for ISR and 
FSR, induce a change in asymmetry through kinemat¬ 
ics. The NLO contributions to the QED interference 
are studied with the latest version of the PHOKHARA 
generator ng, which includes a full matrix element com¬ 
puted at NLO. NLO contributions are found to affect the 
charge asymmetry by 1-2x10“^ in the mass range cov¬ 
ered by the present analysis, where events are generated 


with the highest energy photon in the 20° < 0* < 160° 
reference range [5g . The PHOKHARA results with fully 
implemented NLO corrections are consistent with the 
AfkQed results with independent extra photons, with 
some small discrepancy up to 10“^ for masses larger 
than 4GeV/c^. This shows that the small difference 
between LO and NLO asymmetry originates essentially 
from kinematic effects due to the extra photon. 

The contribution from exchange is investigated 
with the KKMC generator [30], either processed in the 
QED-only configuration, or including the full j + Z^ ex¬ 
change diagrams. As in AfkQed, extra photons are gen¬ 
erated independently in the initial state and final state 
(with PHOTOS). Electroweak (EW) effects are found to 
be at a few 10“^ level, averaging over the full mass range. 

A significant difference of (0.81 ± 0.16) x 10“^ is ob¬ 
served between the asymmetry slopes in KKMC and 
AfkQed, with an asymmetry slope Aq larger (in absolute 
value) in AfkQed than in KKMC. The conclusion holds 
if one considers the forward and backward hemispheres 
separately. As already observed for AfkQed, the asym¬ 
metries expected from KKMC in the two hemispheres 
are consistent with each other. 

Comparison of the asymmetry slope measured in data, 
after acceptance correction, to the full QED-1-.^ ex¬ 
pectation, as implemented in KKMC, confirms that a 
significant difference of 2x10“^ remains in the 1.5- 
4GeV/c^ mass interval, mostly in the forward hemi¬ 
sphere {cos 9* > 0). 
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1.5~4GeV 



VII. RESULTS ON THE CHARGE 
ASYMMETRY IN THE e+e" ^ 77+77-7 
PROCESS 


The charge asymmetry for e+e” —> 77 + 77-7 data be¬ 
fore and after backgroun d su btraction and efficiency cor¬ 
rections is shown in Fig. 21 As for e+e- —>■ /r+zi-y, the 
background for the 77 + 77-7 process is estimated with 
MC, as explained in Sec. IV E The overall efficiency is 
obtained with full simulation of e+e- —>■ 77 + 77-7 events, 
as a function of cos (jf at respective tttt masses, and cor¬ 
rected for data/MC differences in detector response. 

The slopes Aq of charge asymmetry in various TOttw 
intervals are obtained by fitting the corrected charge 
asymmetry distributions to Aq cos (j)* ■ The results for 
the data are shown in Fig. 22 as a function of In 
the p resonance region, the measured asymmetry is neg¬ 
ative, and its magnitude does not exceed ~ IQ-^. A 
clear interference pattern is observed at higher mass. 


A. Comparison and fit to models 


FIG. 20: The difference between the measured asymmetry in 
the p~'y process and the QED prediction, as implemented 
in AfkQed, as a function of cos (j)* for the events in the mass 
interval 1.5-4 GeV/c^, exclnding the J/ip 3.0-3.2 GeV/c^ re¬ 
gion. The result of a linear fit is shown by the solid line. 


4- Conclusion on systematic uncertainties 


In the large number of tests, both experimental and 
theoretical, that have been performed, the antisymmet¬ 
ric part (Aq) of the charge asymmetry is found to be 
remarkably stable. It is immune to all detector ef¬ 
fects taken into account in the simulation, unlike the 
symmetric part (Bq). The simulation properly corrects 
known effects after data/MC adjustment of separate 
sources. The experimental absolute systematic uncer¬ 
tainties on Aq are estimated to be 0.5x10-^ from MC 
reweighting, 0.3x10-^ from data/MC efficiency correc¬ 
tions, O.TxlQ-^ from acceptance edge effects, which sum 
up to 0.9x10-^. In view of the observed differences on 
Aq using AfkQed, PHOKHARA 9.0 (LO and NLO) and 
KKMG (with and without EW corrections), we conser¬ 
vatively set a l.OxlQ-^ systematic uncertainty on the 
theoretical prediction. Adding experimental and theo¬ 
retical uncertainties quadratically a total absolute sys¬ 
tematic uncertainty of 1.4x10-^ is obtained. 

Although we have been unable to find a bias produc¬ 
ing the observed shape as a function of mass of the dif¬ 
ference between the measured Aq and the QED predic¬ 
tions, all data points are within the estimated systematic 
uncertainty, except for 5 out of 14 points near 3GeV/c^ 
that exceed the systematic uncertainty by about 1-2 sta¬ 
tistical standard deviations. 


The magnitude of the charge asymmetry, and its 
variation with mass, measured in the 77 + 77-7 data 
(Eig. 22) is quite different from the prediction of the 
FSR model 1, which treats the pion as a point-like par- 



fit, while the mass nif^ and width F/^ for the /2(1270) 
resonance are fixed to the world averages m- Because 
the measured charge asymmetry loses precision near the 
7777 production threshold and above 1.4GeV/c^, the fit is 
performed between 0.3GeV/c^ and 1.4GeV/c^. The up¬ 
per limit removes the delicate region around 1.5GeV/c^ 
where the pion form factor has a very pronounced dip 
leading to a poor knowledge of the ISR amplitude. 

A distinctive interference pattern is observed in 
Fig. 22 at the location of the /2(1270) resonance. In 
Eq. (12), assuming the dominance of helicity 0 for the /2 


production, the angular dependence of the interference 
term in the tttt c.m. is given by the Legendre polyno¬ 
mial P 2 (cos 0 *), which changes sign at | cos0*| = I/'/i. 
As a consequence, the charge asymmetry is expected to 
follow the same pattern in the vicinity of the /2(1270) 
resonance. 

To check this feature, the charge asymmetries are 
measured separately in the phase space below and above 
I cos 9*\ = 1/ 7 / 3 . The data sample is split accord¬ 
ing to the additional requirement | cos 9*\ < 1/v^, or 
I cos 9*\ > 1/-x/S. To enhance the efficiency in the high 
I cos 9* I region, the event selection is loosened, by re¬ 
moving the p > 1 GeV/c and pion identification require¬ 
ments on the track with lower momentum. To keep 
backgrounds at manageable levels, the higher momen- 
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FIG. 21: Charge asymmetry for e^e —>■ tt+tt 7 data before (•) and after (o) efficiency corrections, in 0.1 GeV/c^ mass intervals 
from 0.3GeV/c^ to 1.8GeV/c^. The line shows the result of the fit to Aocos(f>*. 
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FIG. 22: The fitted slope Aq of the charge asymmetry for 
e^e“ —^ 7r"'"7r“7 at = 10.58 GeV, corrected for efficiency 
within the 20° < 6* < 160° acceptance in the e''"e“ c.m. 



FIG. 23: The charge asymmetries measnred with |cos0*| 
below and above l/\/3. The curves represent the fit results 
(see text). 


turn pion is required to satisfy the tighter identihca- 
tion criteria of a ‘hard tt’ [6]. To further reduce the 
electron contamination, an enhanced Ecai/p < 0.6 se¬ 
lection is applied to the high momentum track, and 
the ionization energy loss in the DCH of the low mo¬ 
mentum track is required to be below the average elec¬ 
tron loss (dS/d^DCH < 650). The corresponding ef¬ 
ficiencies are obtained separately from the full simu¬ 
lation of e“'"e“ —>■ 7r+7r“7 events in the low and high 
I cos 0*1 regions. While the effective |cos0*| range is 
limited to | cos0*| < 0.8 with the standard event selec¬ 
tion, the specific selection applied in the high | cos 9* \ re¬ 
gion allows to extend the asymmetry measurement up to 
I cos 0*1 ~ 0.95. Backgrounds are estimated accordingly 
using the full simulation of relevant processes. After 
background subtraction and overall acceptance correc¬ 
tion, the charge asymmetries obtai ned for | cos 0* | below 
and above l/-\/3 are shown in Fig. 
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Since a large fraction of the events in the standard 
analysis are in the low | cos 0* | region, the charge asym¬ 
metry measured with | cos 0* | < 1/ y/S is quite close 
to the one obtained using the full sample (Fig. 22). 
Although limited by statistics, the charge asymmetry 
measured in the high |cos0*| region presents the op¬ 
posite sign oscillation around the /2(1270) mass, which 
is the expected pattern. The change of sign between 
I cos 0*1 < l/-\/3 and |cos0*| > l/-v/3, and the oppo¬ 
site variation across the resonance provide a solid vali¬ 
dation that the observed charge asymmetry around the 
/2(1270) resonance is indeed due to the interference be- 

+7r“ and 


tween the two amplitudes for e“''e 




7lSR^ 


7FSR/2(1270)(7r+7r ), with /2 in the helicity 0 


state. In the mass range below 1 GeV/c^, the asymme¬ 
try keeps the same (negative) sign in the two | cos 0* | 
regions, as expected for the interference with a scalar 
amplitude flat in | cos0*|. 

The two independent data samples, with | cos 0* | be¬ 
low and above 1 / '/S, are fitted separately to the model, 
and the fitted cq and C 2 amplitudes are obtained in both 
cases. Since the pure ISR AfkQed MC used to compute 
the efficiencies is not expected to properly correct for 
unmeasured regions of | cos 0* |, the fit of the data above 
l/-\/3 is performed in the effective range of non-null ef¬ 
ficiency where asymmetries are measurable. 

The model describes the data well, and the two sets 
of htted values of cq and C 2 are consistent in sign and 
magnitude and can be averaged, yielding cn = —1.27 ± 
0.20 and Ca = 5.4 ± 1.6. 


B. Monte Carlo reweighting and final results 


Since a significant asymmetry is observed in the data 
in contrast with AfkQed, which does not include LO 
FSR in the pion channel, the overall efficiencies e± com¬ 
puted with MC and used to measure the asymmetry in 
data are biased through the cross effect be tween t he ac¬ 
ceptance and the physical asymmetry (Sec V A 1). This 
situation calls for an iterative procedure to introduce the 
observed interference effect into the MC. 


To implement this procedure, new MC samples of 
reweighted events are produced, in which the weights 
are computed event by event as the full cross section 
value including LO FSR divided by the ISR-only cross 
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FIG. 24: The charge asymmetry slopes for e^e“ —>■ 7r+7r“7 using MC samples with and without reweighting, and fit to the 
model after reweighting (see text); (top left) for | cos0*| < l/v^; (top right) for | cos6^*1 > l/v^; (bottom) for the full | cos0*| 
range; the blue band represents the model-2 prediction using the average Co ,2 values after reweighting. The light-blue part 
corresponds to the extrapolation of the model beyond the fitted range. 


section, for the values of cos 6*, cos 9* and cost/)* 
for the event. The differential cross sections are given by 
the model used to fit the data (Eqs. (8p0)). The FSR 
model is made quantitative by using the fitted values for 
Co, 2 - The studies have been performed separately for the 
regions below and above l/-\/3. 


The fitted values of Cq and C 2 are stable after two iter¬ 
ations. A third iteration is performed in order to check 
the stability of the results. The difference between the 
last two iterations is taken as a systematic uncertainty. 

The final Aq values are given in Fig.[^ together with 


the FSR model prediction using the fitted co ,2 values de¬ 
termined by the iterative process. The extrapolation of 
the model beyond the fit region 0.3-1.4 GeV/c^ is shown 
(light-blue band). Although the statistical uncertainty 
of the data is large, there is evidence that the model 
becomes inadequate above 1.8GeV/c^. This is not sur¬ 
prising since a constant S-wave amplitude and the /2 
resonance are likely to be insufficient to describe this 
region, where many high mass resonances contribute to 
the 7r+7r“ final state. However, based on the change 
of asymmetry at the 10“^ level induced at lower mass. 
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the effect of inadequate reweighting in the last few mass 
bins is expected to be much smaller than the statistical 
uncertainty. 

Since the two independent sets of parameters agree 
within their uncertainties, they can be combined and 
the weighted average of the fitted values represents the 
best information which can be obtained from this in¬ 
terference analysis (Table H- An alternative is to fit 
the overall sample obtained with the standard selec¬ 
tion, using efficiencies calculated with the reweighted 
MC. In this case, as the data spans over the sign change 
at I cos 0 *\ = 1/-s/S, the measured asymmetry is much 
reduced. Therefore the combined result from the two 
complementary ranges is more sensitive, and moreover 
provides a clear confirmation of the helicity 0 /2(1270) 
contribution. 


C. Systematic uncertainties 

The difference AAq between the results with and 
without background subtraction is found to be less 
than 10“^, except near threshold (1x10“^) and above 
l.lGeV/c^ (1x10“^ at 1.25GeV/c^). Except in the dip 
region at 1.55GeV/c^ where the ISR cross section has 
a sharp minimum and statistical uncertainties are very 
large, the background level has been checked [S] with a 
precision of 20% in the worst cases. The systematic un¬ 
certainty due to background subtraction is consequently 
estimated to be less than 2x10“^ from 0.4 to 1.1 GeV/c^, 
increasing above (2x10“^ at 1.25GeV/c^). 

As observed for e+e“ —>■ the selection re¬ 

quirements for trigger, tracking, and tt-ID have charge- 
asymmetric efficiencies for the e^e“ —>■ 7r+7r“7 process. 
The corrections for the difference between data and MC 
on the efficiencies are included in the overall acceptance. 
The difference between the charge asymmetry results 
with and without the data-MC corrections is smaller 
than 10“^ except in the dip region (0.5x10“^). Since 
the corrections are determined with data with a pre¬ 
cision of 10%, the resulting systematic uncertainty is 
negligible. 

As done for /r+^“ 7 , the asymmetry A{cos(j)*) in 
7 r+ 7 r “7 is alternatively fitted to AqCOS^* -I- Rq since a 
bias on Aq and Bq values inconsistent with zero might 
disclose an incorrect efficiency determination, or an in¬ 
correct background subtraction. As shown in Fig. |25[ 
the fitted slopes Aq deviate from the final values by 
less than 10“^, except in the background dominated dip 
region of the cross section (m-TrTr — 1.5GeV/c^), where 
the deviation AAq reaches 3x10“^. In the mT^TT region 
where the fit of the theoretical model is performed (0.3- 
1.4GeV/c^), the slope Aq does not change by more than 
5x10“^, except in the last (1.3-1.4GeV/c^) bin, where 
the deviation is 2x 10“^. This is consistent with the esti¬ 
mated background contribution to the systematic error. 
The fitted Bq values are within 2.5o’ from zero over the 
full mass range. The average Bq for m,r 7 r < 1 GeV/c^ is 


(0.41 ±0.16) X 10-2. 

Different interaction rates in the detector material for 
positive and negative pions induce a charge asymme¬ 
try. Although such an effect is included in the simu¬ 
lation of the detector response based on GEANT4, its 
description and the corresponding track loss are known 
to be somewhat imperfect. Independent studies have 
shown that data/MC discrepancies occur at the 10% 
level for both tt"*' and tt”, in opposite directions. A 
residual charge asymmetry is thus expected after ap¬ 
plying the MC corrections. The effect of imperfect 
simulation of nuclear interactions is investigated using 
the large sample of Tr+Tr-y events produced by AfkQed 
at the generator level. A weight is assigned to each 
track according to its momentum and its path length 
through detector material as a function of the polar an¬ 
gle, using a ±10% relative change in the respective 7r+ 
and Tr~ interaction rates. The charge asymmetry ob¬ 
tained after the interaction reweighting is subjected to 
the two-parameter linear fit in the integrated mass range 
from 0.4 to 1.2 GeV/c^. The slope Aq changes by only 
AAg = (—0.006 ± 0.024) x lO-^, which confirms the ro¬ 
bustness of the Aq observable. The charge asymmetry it¬ 
self is however modified as the fitted Bq value is found to 
be displaced significantly, ABq = (0.240±0.016) x lO-^, 
in good agreement with the observed Bq value in data 
in the same mass range. Imperfect simulation of nuclear 
interactions thus provides a plausible explanation of the 
small Bq values found in the analysis, while leaving the 
Aq measurement unaffected. 

Summing up all sources, including the estimated 
cross-effect between acceptance and physical asymme¬ 
try, the absolute systematic uncertainty on Aq is esti¬ 
mated to be less than 0.17% in the /2(1270) region and 
less than 0.1% elsewhere. 


D. Searching for an / 2 ( 1270 ) signal in the 7r“*“7r 7 
cross section 

Given the sizeable amplitude C 2 of the D-wave con¬ 
tribution to the ISR-FSR interference obtained from 
the charge asymmetry measurement, direct evidence of 
/2(1270) production is searched for in the cross section 
measurement. While the latter is overwhelmingly dom¬ 
inated by the ISR production of the p resonance, the 
rapid fall off of the pion form factor in the vicinity 
of the /2(1270), and the distinct angular distribution 
P2 (cos 0*) of the D-wave in the tttt system, are assets 
used in the direct search. 

Since the P 2 {cos 6 *) distribution exhibits a peak at 
I cos 0*1 near unity, in contrast with the sin^ 0* depen¬ 
dence of the ISR cross section, the search is performed 
in the very high range | cos 0* | >0.85. Because the stan¬ 
dard event selection depopulates that region completely, 
due to the momenta of both tracks being required to be 
larger than 1 GeV/c, the direct search uses the specific 
selection designed for the charge asymmetry measure- 
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TABLE II: The parameters obtained from the fit of the charge asymmetry for e^e“ —>■ 7r"'"7r“7 at = 10.58 GeV after three 
iterations, with 20° < 9* < 160° in the e''"e“ c.m., in the |cos6*| below and above l/%/3 regions, and the weighted average, 
where the errors are statistical. The results of the direct fit over the full range are given in the last column. 


parameter 

1 cos 0*1 < 1/-\/3 

1 cos 0*1 > l/v^ 

average 

all 1 cos 0* 1 

Co 

-0.84 ±0.24 

-1.13 ±0.35 

- 0.93 ± 0.20 

-0.87 ±0.20 

C2 

3.82 ± 1.81 

6.33 ±3.03 

4.48 ± 1.56 

3.41 ±4.25 


< 0.03 

0.02 


0.01 


0.5 

1 

1.5 2 

0.5 

1 

1.5 2 



(GeV/c^) 



(GeV/c) 




FIG. 25: The change of slope AAq (left) and constant term Bo (right) in the fit of charge asymmetry to Aq coscp* + Bo, as a 
function of tttt mass for e''"e“ —>■ 7r+7r“7 . 


ment in the high |cos0*| region (Sec. VIIA I, with an 
even tighter dE/dxDCH < 550 requirement. 

The TTTT mass spectrum of the reconstructed events 
in this specific analysis is displayed in Fig. with the 
largest expected background from /x/ry events to be sub¬ 
tracted. The resulting spectrum is dominated by the 
ISR production, which is also subtracted. The remain¬ 
ing spectrum shown in Fig. does not present any sig¬ 
nificant excess at the /2(1270) mass or elsewhere, ex¬ 
cept for a slow rise above 2 GeV/c^ (not shown) that 
originates from a residual eey background. The mass 
spectrum is fitted between 0.95 and 1.95GeV/c^ to a 
constant and a Breit-Wigner lineshape with the world 
average /2(1270) mass and width [31] and a free-floating 
amplitude. The fitted number of /2 events in the mass 
interval at the peak is found to be 4.7 ± 4.2, to be com¬ 
pared to 26.7±1.1 ISR events in the same interval. After 
correction for the loss of efficiency near | cos 0*1 = 1 ob¬ 
tained from MC for ISR and /2 candidates, the /2(1270) 
fraction |/ 2 p/(|ISRp -b I/ 2 P) in the / 2 -enhanced range 
0.8 < I cos 0*1 < 1 is measured to be 0.22 ± 0.15. This 
corresponds to a |c 2 | value equal to 4.6 ± 2.2. 

The three independent determinations of |c 2 | (the 
interference fits in two cos 9* regions and the direct 
72 ( 1270 ) search in the cross section) yield consistent re¬ 


sults. Since a positive sign is clearly indicated by the 
interference analysis, the value from the direct search is 
also taken to be positive. The three independent values 
can be combined with the result C 2 = 4.5 ± 1.3, estab¬ 
lishing LO FSR production of the /2(1270) resonance at 
the 3.6 ct level. The corresponding production cross sec- 
tion is (377 i 8) fb. The results are displayed in Fig. 

The size of C 2 is about a factor of two larger than 
value predicted by Ghernyak |32| with a QCD model giv¬ 
ing |c 2 ^| = 2.2. However the difference only amounts to 
1 . 8 ( 7 , not including the unknown theoretical uncertainty. 
The sign is not provided in Ghernyak’s prediction. 
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E. Consequences for the cross section 
measurement by BABAR for e“*“e“ -A- tt+tt" and 
contribution to the anomalous magnetic moment of 
the muon 

In the measurement of the e"'"e“ -A cross 

section by the BaBAR collaboration [Bj using the ISR 
method, the lowest-order FSR contribution was argued 
to be negligible, based on theoretical estimates. The 
primary result of the present interference analysis is to 
determine the actual size of the |A4 fsrP cross section. 
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FIG. 26: The tttt mass spectrum for | cos6*| > 0.85 for the 
analysis extended to low momentum (points), the expected 
background from misidentified events (red histogram), 
and the predicted ISR spectrum from the standard cross sec¬ 
tion analysis (blue histogram). 



FIG. 27: The tttt mass spectrum for |cos0*| > 0.85 for the 
analysis extended to low momentum (points) after subtrac¬ 
tion of the fj.jj.'y background and the ISR contribution. The 
curve is the result of a 6t to a constant term and a simple 
Breit-Wigner shape for the /2 resonance. 


misinterpreted as ISR, and its contribution to the total 
cross section (|A^isrP + |A^fsrP)- 

Using the FSR model 2, which describes well the mea¬ 
sured charge asymmetry in the [0.3-1.4] GeV/c^ range, 
the FSR cross section calculated through Eq. (§ with 
the fitted co ,2 parameters, is extrapolated to higher 
masses. The resulting FSR fraction in the BABAR cross 
section is given in Fig. as a function of m^ 7 r- As 
expected the FSR fraction is negligible in the p re¬ 
gion, but increases significantly above 1 GeV/c^ due to 
the /2(1270) contribution and the rapid fall-off of the 
pion form factor. In fact the FSR ‘background’ exceeds 
the estimated systemati c un certainty quoted in Ref. [B] 
(green histogram in Fig. 291 for mass above 1.2GeV/c^, 
while remaining close to the total uncertainty (black his¬ 
togram) . The FSR contribution is found to be dominant 
around 1.5GeV/c^, in the region where the ISR cross 
section displays a deep dip and is consistent with zero 
within the large errors. 


The contribution to the muon magnetic anomaly, 
= (5/i ~ 2)/2, from hadronic vacuum polarization 
involves a dispersion integral over the cross section 
e“'"e“ —>■ hadrons weighted by a known kernel (Ref. |33j 
and references therein). The integral is dominated by 
the 7r'''7r“ channel and its most precise determination to 
date is from BABAR using the ISR method [B] with the 
value = (514.09 ± 2.22,tat ± 3.1Uyst) x 10-1° 

when integrating from threshold to 1.8 GeV. This value 
is derived under the assumption that the cross section for 
e“'"e“ —>■ 7r'''7r“7(7) has a negligible contribution from 
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FIG. 28: The results obtained for the /2(1270) amplitude C 2 
in the process e''"e“ —>■ 7r'*’7r~7 from the interference analysis 
and the direct /2(1270) search in the cross section. The labels 
‘low’ and ‘high’ refer to the determination in the two angular 
ranges |cosS*| < l/\/3 and l/\/3 < |cos6*|. For the direct 
search the positive solution is chosen. The combined value 
for the three independent analyses is given by the vertical 
band. 
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FIG. 29: The FSR fraction in the BABAR measurement [B] of 

the 7r7r7 cross section, defined as the ratio -nr;—- rr 

obtained in this analysis using FSR model 2 with co ,2 pa¬ 
rameters fitted to data (blue band). The light-blue part cor¬ 
responds to the extrapolation of the model beyond the fitted 
range. The FSR fraction is compared to the systematic er¬ 
ror of the BABAR cross section measurement (green dashed 
histogram) and its total error (black histogram). 


LO FSR. 

The present measurement of the charge asymmetry 
allows one to validate this assumption in a quantitative 
way. Using the FSR fraction shown in Fig. the con¬ 
tribution to Qfj, from the LO FSR falsely attributed to 
the ISR cross section is found to be, in the same energy 
range up to 1.8 GeV 

Aal'^iFSR) = (0.26 ± 0.12) x 10”^°. (17) 

This reduces the value of by (5.1±2.3)x 10“"* relative 
to the BaBAR determination. The correction is small 
compared to the total BaBAR relative uncertainty of 7.4x 
10“^, which justifies its earlier neglect. 


VIII. CONCLUSIONS 

The radiative process e"'"e“ —)■ Xj, where X = 
and TT~^TT~ are considered in this analysis, involves con¬ 
tributions from both LO ISR and FSR. Because charge 
parities of the final state pair are opposite for ISR and 
FSR, the interference between ISR and FSR changes 
sign with the charge interchange of the two muons (pi- 
ons). As a consequence, investigation of the charge 
asymmetry of the process gives a way to study the in¬ 
terference between ISR and FSR, which is sensitive to 
the relative contribution of LO FSR. 


From QED for and from FSR models for 

7 r+ 7 r“q, we find that the charge asymmetry A has a 
strong dependence on the angle (j)* between the 
{Tr'^'!T~) plane and the e“'"e “7 plane in the e“'"e“ c.m. 
system, which can be simply represented by a linear 
function A = Aq cos (j)*. The slope Aq quantifies the 
magnitude of the interference between ISR and FSR. 

The acceptance effects on the measured charge asym¬ 
metry are studied with the full simulation of e“'"e“ —)■ 
rArT^ and e“''e“ —> 7r+7r“7 events. We find that the 
detector and event selection, including trigger, tracking, 
PID, and kinematic fitting, induce nonlinear patterns 
on the cos (jf dependence of the charge asymmetry, but 
have a small impact on the determined slope Aq. Kine¬ 
matic acceptance — namely the angular acceptance, and 
energy or momentum requirements on the final state 
particles — changes the slope of the observed charge 
asymmetry significantly, although the kinematic require¬ 
ments are charge-symmetric. This is due to a cross ef¬ 
fect between acceptance and true interference that pro¬ 
duces a bias in the measured asymmetry if the physical 
asymmetry differs between data and MC. This bias is 
corrected through an iterative procedure in the Tr+Tr^y 
analysis, as in that case the charge asymmetry is null in 
the generator. 

After background subtraction and correction for the 
overall acceptance, which are obtained from the full sim¬ 
ulation with corrections for data/MC differences, we 
measured the slope Aq of the charge asymmetry as a 
function of (m,r 7 r)- The QED test, namely the com¬ 
parison between the charge asymmetry measured in the 
rA rT ^ data and predicted by the simulation, in which 
the LO ISR-FSR interference is implemented, shows an 
overall good consistency. However, some absolute de¬ 
viation amounting to AAq = ~ 0.03 in 

the 3GeV/c^ region is observed and cannot be fully ex¬ 
plained by known systematic effects, either in the data 
or in the MC generators, which are estimated to be less 
than 0.014. 

The measured slope Ag of charge asymmetry in the 
e+e” —?► 7r+7r“7 data is about —1% and flat around 
the p mass. Outside of the p peak, the data exhibits 
the pattern expected from the interference between 
e+e” —>■ 7isR7r^7r“ and e+e“ —)■ 7FSR/2(1270)(7r+7r“). 
The data shows a good consistency with the predictions 
of a model of ESR from quarks with contributions of 
a scalar widespread mass distribution and the /2(1270) 
tensor resonance. In the p region the results are not 
consistent with a model based on ESR from point-like 
pions (scalar QED), in contrast with the observations at 
low energies [4]. 

These results are first measurements of the charge 
asymmetry in the e“'"e“ —>■ pApT^ process, and for 
e+e“ —)■ 7r+7r“7 at high energy {^/s ~ 10.58GeV). The 
FSR contribution to e“'"e“ —>■ 7r“'"7r“7 derived from this 
analysis is small and this confirms that it is negligible in 
the measurement of the cross section obtained by BaBAR 
assuming pure ISR [B]. Accordingly this FSR bias trans- 










32 


lates into a correction to the muon magnetic anomaly of 
only (0.51 ± 0.23) per mille of the tttt hadronic vacuum 
polarization determined from BaBAR data, small com¬ 
pared to the total quoted uncertainty of 7.4 per mille. 
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